Short Note 



Typeset with jpsj2.cls <ver.l.l.l> 



Decoherence of interacting electrons in disordered conductors: on the relation 
between influence functional and diagrammatic approaches 

Jan von Delft* 

Lehrstuhl fur Theoretische Festkorperphysik, CeNS und Sektion Physik, LMU Miinchen, Theresienstr. 37, D-80333 Miinchen 

(Received August 7, 2002) 

We establish a connection between the influence functional approach of Golubev and Zaikin 
(GZ) and Keldysh diagrammatic perturbation theory for calculating the decoherence time t v of 
interacting electrons in disordered metals; we show how the standard diagrams for the Cooperon 
self energy can be recovered from GZ's influence functional e ~(' s R +s '\ This allows us to shed 
light on GZ's claim that Sr is irrelevant for decoherence: Sr generates as many important self 
energy diagrams as Si; GZ's neglect of Sr is permissible only at high temperatures (T > h/r e \). 



KEYWORDS: Decoherence, influence functional, diagrammatic perturbation theory, interactions, disorder 



1. Introduction 

A few years ago, Golubev and Zaikin (GZ) developed 
an influence functional approach for describing interact- 
ing fermions in a disordered conductor. Their key idea 

' was as follows: to understand how the diffusive behav- 
ior of a given electron is affected by its interactions with 

' other electrons in the system, which constitute its ef- 

, fective environment, the latter should be integrated out, 
leading to an influence functional (denoted by e~ < ^' lS^+SI ' , 

, below) in the path integral J T>(RP) describing its dy- 

■ namics. To derive the effective action iSr + Sj, GZ 
devised a strategy which, when implemented with suf- 

■ ficient care, properly incorporates the Pauli principle - 
, this is essential, since both the particle and its environ- 

■ ment originate from the same system of indistinghuish- 
able fermions, a feature which makes the present problem 

' interesting and sets it apart from all other applications 
of the influence functional strategy that we are aware of. 
GZ used their new approach to calculate the electron 

, decoherence time t v (T) , as extracted from the magneto- 
conductance in the weak localization regime, and found 
it to be finite at zero temperature: 1-4 ) t v (T — ► 0) = r°, 
in apparent agreement with some experiments. 5 ^ How- 
ever, this result contradicts the standard view, based 
on the work of Altshuler, Aronov and Khmelnitskii 6 ) 
(AAK), that t v (T — > 0) = oo, and hence elicited a 
considerable and ongoing controversy, 7 '' with pertinent 
critique coming, in particular, from Ref. 8,9 -' 

The fact that GZ's final results for t v (T) are contro- 
versial, however, does not imply that their influence func- 
tional approach, as such, is fundamentally flawed. To the 
contrary, having repeated their calculations in detail, we 
have come to the conclusion that their strategy is sound 
in principle and that an influence functional of the form 
e -(iS R +Si) which they found can indeed be derived with- 
out making non-standard approximations. In fact, it can 
be shown, and this is our main result, that the standard 
Keldysh diagrammatic expressions for the self energy of 
the Cooperon can be obtained from iSr + Si. However, 
when applying this influence functional to the problem 
of decoherence, GZ make a semiclassically-motivated ap- 



proximation according to which the the effects of Sr 
for decoherence can be neglected. We recount a sim- 
ple back-of-the-envelope argument, 10 - ) due to F. Mar- 
quardt, to suggest that neglecting Sr is permissible only 
for T > h/r e \, where r e i is the elastic mean free time. 
We shall also show that in diagrammatic language, ne- 
glecting Sr corresponds to neglecting several important 
diagrams contributing to the Cooperon self energy, as 
first pointed out in Ref. 9 -* 

Equation numbers from GZ's papers, 1-4 ^ will be pref- 
aced, when cited below, by I, II, III or IV, respectively. 

2. The model 

We consider a disordered system of interacting 
fermions, with Hamiltonian H = Hq + Hi, where 

H Q = J dxft(x)ho(x)$(x) , (1) 

R\ = ^-Jd Xl dx 2 :ft{x l )j>{x 1 ):Vi? :ft(x 2 )i>(x 2 ): 

Here J dx = J dr, and tp(x) = ip(r, a) is the electron 
field operator for creating a spin-tr electron at position 
r, with the following expansion in terms of the exact 
eigenfunctions ^(a;) of h {x) = ^-V^ + V imp (r) - fx: 

i){x) =^ A (a;)c A , [ftoW-WxW = 0. (2) 

A 

The interaction potential = V^flri — r 2 \) acts be- 
tween the normal-ordered densities at r\ and r 2 . 

3. Influence functional for interacting electrons 

GZ proposed a strategy (whose steps and approxima- 
tions are recapitulated in the Appendix) , that allows the 
DC conductivity, ctdc, to be expressed [see Eq. (A-3)] in 
terms of path integrals of the following general form [cf. 
(11.53), (IV.31)]: 

C 12 /,3y = I F B V(RP)e-^ §R+Sl ^' t2 ^ h . (3) 

J2' p ,2b 
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Fig. 1. A pair of paths contributing to the weak localization cor- 
rection, illustrating Eq. (A-3). The interaction lines, gener- 
ated by expanding the influence functional e~ ( lS R+ s i) as in 
Eq. (A-14), are labelled according to Eq. (A-13). 



interaction propatators (R — C R , I — i\C K ). 

Via the influence functional, Eqs. (3) to (6) concisely 
incorporate the effects of interactions into the path in- 
tegral approach. Si describes the classical part of the 
effective environment, and corresponds to the contribu- 
tion calculated by AAK. 6 ' With Sr, GZ succeeded to 
additionally also include the quantum part of the envi- 
ronment, and in particular, via the occurence of the den- 
sity matrix (R, P) in Eq. (6a) , to properly account for 
the Pauli principle. Note, though, that Eqs. (3) to (6) all 
refer to a given impurity configuration; impurity averag- 
ing still has to be performed, and in the path integral 
formalism it is by no means easy to do this properly. 



The symbol J T>(RP) is a shorthand for the following 
coordinate-momementa double path integral, 



I 



2' F ,2 B 

H B (ti)=»V 



R F (ti)=n 

R F (t 2 )=r 2 , 



V{RP) = / VR F (t 3 ) I VP F {t 3) 



x / 'VR B {t 3 ) I VP B (t 3 )e^ S o- S o](^)^, 

JR B (t 2 )=r 2 J 

(4) 

which, when taken by itself, gives the amplitude for a 
free electron to propagate from r 2 > at time t 2 to n at 
ti, times the amplitude for a free electron to propagate 
from rv at time t\ to r 2 at t 2 [corresponding to the loop 
parts of the paths in Fig. 1] , in the absence of interactions 
with other electrons. We shall call these the forward and 
backward paths, respectively, and label them by an index 
a = F,B. The corresponding free actions S§ = 5 are 
given in Eq. (A-4). The weak localization correction to 
the conductivity, o~^q , arises from contributions to odc 
for which the coordinates T\, r[, r 2 and r' 2 all lie close 
together. We henceforth consider only this case. Then 
£-12'. 21' i s j us t the Cooperon propagator, dominated by 
contributions from those classical paths for which path 
B is the time-reversed version of path F . The effect of 
the other electrons on this propagation is encoded in the 
influence functional e~^ lSR+Sl ^ occuring in Eq. (3). The 
effective action iSn + Si turns out to have the form 



SR/i{ti,t 2 )= ^2 [ dt 3 [ dt ± L zll a , 



(5) 



where the L 3a4 , are functions of the coordinates and 
momenta R a (t 3 ), P a (t 3 ) and R a ' (i 4 ), P a ' (t 4 ) that occur 
in the path integral [cf. (11.54,55)]: 

L$ a4a , =s a ,R(t 3 - t 4 , R a (t 3 ) - R a ' (t 4 )) 

x ±[l-2p a (R a '(t 4 ),P a '(t 4 )] , (6a) 

L Li a , =*aSa>l(t 3 - U,R a (t 3 ) - R a \t 4 )) , (6b) 

Here s a stands for s F / B = ±1, p^(R 7 P) is the single- 
particle density matrix in a mixed position-momentum 
represenation [cf. Eq. (A-6)], while R(t,R) and I(t,R) 
are real functions [given by Eq. (A-7)] that are, respec- 
tively, proportional to retarded and Keldysh parts of the 



4. GZ's strategy for determining t v 

To calculate the decoherence time t v , GZ argue as fol- 
lows: the effective action in Eq. (3) in general causes 
the Cooperon to decay with increasing time, say as 
Ci2',2i' ~ C^ ^e-W* 1 -* 2 ) [cf. (11.66), (IV.10)], where 
fd(t) is an increasing function of time; t v is the time 
scale characterizing this decay, set by fd{T v ) — 1- To ob- 
tain the function fd, GZ evoke a standard semiclassical 
argument: since the path integral is dominated by the 
saddle point paths of the free action S$, i.e. by the set 
of classical, time-reversed diffusive paths, they take fd 
to be (i) the disorder average 1 ^ ( )di s of the sum over 
all classical paths ( ) c ] of the effective action evaluated 
along such a path (iS R l + Sf), but (ii) without including 
any non-classical paths: 



fd(h - h) 



Sf)(h,t 2 )) c 



(7) 



[cf. (III. 22), (IV. 11)]. Moreover, in the spirit of semiclas- 
sical approximations, they (iii) take the limit 12 ' h — > 
in the single-particle density matrix occuring in Eq. (6a) 
for Sr, i.e., they replace Pq(R,P) of Eq. (A-6) by 
the "occupation number" n (ho{R, P)) , where no(£) = 
l/[e«/ T + 1] is the Fermi function [cf. (11.43), (11.68)]. 

Within the approximations (i), (ii) and (iii), GZ find 
that S R l ~ for any given pair of classical, time-reversed 
paths, and hence conclude that Sr is "irrelevant" for 
decoherence, which is thus determined by Si alone [cf. 
discussion before (III. 22), or after (IV. 31).] They thus 
calculate fd purely from Sf [see, e.g., (IV. 12)], and 
from fd{T v ) — 1 find, in 1 dimension for example, [cf. 
(11.76,77)] 



1 

Tin 



duj coth(w/2T) 




(8) 
(9) 



Note that the frequency integral has an ultraviolet diver- 
gence at large u>, and hence has to be cut off by hand. 13 ' 
GZ cut it off at w^ax = V T cb the inverse elastic mean 
free scattering time, arguing 13 ' that at higher lo (smaller 
times) the "approximation of electron diffusion becomes 
incorrect" [cf. paragraph before (11.76)]. This leads to a 
finite decoherence time at zero temperature in Eq. (9), 
^/t<p{T — > 0) 7^ 0. In contrast, according to the phi- 
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losophy of AAK, 6 ) one should take w** K = T/h, since 
frequencies larger than T/h would correspond to virtual 
excitations of the environment, which are believed not to 
contribute to decoherence. This would yield the standard 
result 1/t v (T — > 0) = 0. Thus, the controversy centers 
on the question whether the higher frequency modes do 
contribute to low temperature decoherence or not. 

5. On the importance of S R for decoherence 

Although we agree with GZ's influence functional (3) 
to (6), we disagree with GZ's central conclusion that Sr 
is irrelevant for decoherence. Firstly, influence function- 
al have the general feature that deducing decoherence 
properties from Sj alone is reliable only for highly ex- 
cited states: as emphasized in a very insightful recent 
paper by F. Marquardt, 14 ) neglecting Sr would amount 
to neglecting dissipative effects by which the quantum 
system can give back to the environment some of the 
energy which it picks up from the classical part of the 
environment; this neglect of dissipation would thus cause 
the system to heat up beyond what is allowed by detailed 
balance (an effect which can be ignored for decoherence 
only if the system is already highly excited). 

Secondly, and more specifically: even if one puts aside 
all reservations 11 ' 12 ^ about approximations (i) and (iii) 
[although these reservations are serious], and accepts 
GZ's result that Sr ~ for classical paths, a question 
remains about approximation (ii): what about fluctu- 
ations, i.e. contributions from paths in the path inte- 
gral that deviate slightly from the classical paths"? If 
the "leading" contribution vanishes, the "sublcading" 
one becomes of interest! Indeed, it is natural to ex- 
pect that such close-to-classical paths can produce a sig- 
nificant contribution missed by GZ's semiclassics, say 
SSr =/= 0, due to the following intrinsic problem (IP): 
Sr depends on the density matrix po(R,P), which at 
low termperatures is a very sharp function of R and P, 
and hence very sensitive to fluctuations of ( or approx- 
imations inside 12 ) ) its arguments. This intrinsic prob- 
lem sets the present problem apart from other, exactly 
solvable models such as the Caldeira-Leggett model, and 
renders meaningless any attempts 3 ' to justify neglecting 
Sr here by comparisons to such other models. 

Let us estimate 10 - 1 the importance of these fluctua- 
tions, by estimating SSr for a segment of a diffusive path 
in which an electron travels ballistically with velocity vf 
(and energy ~ fi) for a time r e i between two scatter- 
ing events at two impurities. Now, the close-to-classical 
paths that still contribute significantly to the path inte- 
gral are those for which [So — Sq 1 ) < h. Considering only 
the first term in Eq. (A-4) for So, we thus use 



2'f R 4 f 




3f R 1f 






3 b A 1 b 


2b A 4 b 



ft ~ (So - S^ 1 ) ~ [ * dtSP-R~ r e] 5Pv F , 
Jo 



(10) 



or SP ~ ft/T e iUF (which is just the uncertainty princi- 
ple), to estimate the typical magnitude of momentum 
fluctuations around Pp. It follows that 



Sn (h {R, P)) 



5h (R, P) _ P F SP _ _h_ 
T ~~ mT ~ ^T' 




(b) (22*) 
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Fig. 2. First order contributions to the irreducible self energy of 
the Cooperon, illustrating (a) Eq. (A- 15), and (b,c) Eqs. (A- 16). 
(d): typical path contributing to S^ F ; GZ neglect T, R and hence 
such paths, as first pointed out in Ref. 8 ' 



Now, GZ neglect SSr relative to Sf; since this requires 
Sn a < 1, it can be reliable only for T > h/r e i. For smaller 
T, the contributions of SSr should become important. 15 ) 
Having realized that SSr is important at low temper- 
atures, the natural next question is: can SSr cancel the 
ultraviolet divergence arising from Sf? We believe it 
does, but showing this will require a more accurate cal- 
culation than GZ's, that does not resort to semiclassical 
arguments and avoids approximations (i) to (iii). 

6. Obtaining diagrams from influence functional 

As a first step in that direction, we have explored the 
connection between GZ's influence functional and stan- 
dard diagrammatic perturbation theory. The connection 
turns out to be remarkably simple: upon performing the 
momentum integrals in J V(RP) and expanding the re- 
sulting influence functional (some details are given in 
App. A. 2), one generates a Dyson-like equation for the 
Cooperon [Eq. (A-14)], with a self energy whose lowest 
order 16 -* irreducible diagrams [given by Eq. (AT6)] are 
depicted diagrammatically in Fig. 2. Remarkably, the 
resulting diagrams coincide precisely with those obtained 
by standard Keldysh diagrammatic perturbation theory, 
as depicted, e.g., in Fig. 2 of Rcf. 9 ) (There, impurity lines 
needed for impurity averaging are also depicted; in our 
Fig. 2, they are suppressed). This fact, which is our main 
new result, is a strong indication that the expressions of 
Eqs. (3) to (6) for the influence functional are sound, 
and the approximations made during its derivation rea- 
sonable [steps (5) and (6) in App. A.l, and Ref. 17 )]. 

Moreover, this fact also allows us to shed new light on 
the roles of Sr and Si in the effective action. Sj gives 
rise to the terms Sr, which contain a factor coth(ftw/2T) 
[arising from the factor I(k 7 uj) 7 cf. Eqs. (A-7b) and 
(AT3)], and Sr gives rise to the terms Ejj, which con- 
tain a factor tanh[(eA — w )/27 1 ] [arising from the factor 
(5 - 2p°) in Eq. (A- 13)]. Now, the sum of all the self- 
energy diagrams for iT,r + E# in Fig. 2 has been eval- 
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uated in Ref., 9 ) in the energy-momentum representation 
and for those choices of energy and momentum argu- 
ments that are relevant for determining the Cooperon 
lifetime; the result [Eqs. (3,4) of Rcf. 9 )] was found to be 
not ultraviolet divergent, since the coth and tanh func- 
tions always occur in the combination 

/j 
— {tanh[(e A - w)/2T] + coth(w/2T)} . . . , (12) 

so that the frequency integral is cut off at w max ~ T j h, as 
anticipated by AAK. 6 ) We expect this cancellation of UV 
divergencies from iSr and Si to occur not only in first 
order perturbation theory, but in every order, 22 ' since 
the structure, e~( lSR+Sl \ of the influence functional is 
such that the self energy contribution iT,R and S/ always 
occur in the combination iY<R + E/. 16 -* 

Conversely, we conclude that GZ's neglect of SSr cor- 
responds to neglecting the contribution of Sr to the 
Cooperon self energy, i.e., GZ neglect all the diagrams 
of Fig 2(b), as first pointed out in Rcf. 9 ) As argued 
above, we believe that this is allowed only if T > h/r e \. 

7. Conclusion 

Our analysis can be summarized as follows: GZ's influ- 
ence functional stratey is sound in principle; when imple- 
mented with sufficient care, it properly incorporates, via 
Sr, the Pauli principle. However, GZ neglect the latter 
by neglecting SSr during the semiclassical calculation of 
t^. This can work only for large temperatures. 

A complete, first-principle evaluation of t v would be 
obtained if one sums up the Dyson- like equation (A- 14) 
for the Cooperon in the presence of disorder, either dia- 
grammatically or by using path integral techniques, but 
without neglecting SSr. To the best of our knowledge, 
this program has not yet been fully carried out to the 
end. Thus, at least in the eyes of the present author, 
the fate of t v as T — + for the present model has not 
yet been worked out in complete, conclusive detail. At 
present we see no reason to believe, though, that the 
result will disagree with the conclusions of AAK. 

Acknowledgements 

I thank D. Golubev and A. Zaikin, and, in equal mea- 
sure, I. Aleiner, B. Altshulcr, M. Vavilov, I. Gornyi and 
F. Marquardt, for numerous patient and constructive dis- 
cussions, which taught me many details and subtleties of 
the influence functional and diagrammatic approaches, 
and without which I would never have been able to reach 
the conclusions presented above. Moreover, I thank V. 
Ambegaokar for raising and sustaining my interest in this 
problem, and acknowledge illuminating discussions with 
J. Imry, J. Kroha, S. Mirlin, H. Monien, A. Rosch, I. 
Smolyarenko, G. Schon, P. Wolfle and A. Zawadowski. 

Note added in proof: — Very recently, GZ defended 
their theory yet again 23 ) in cond-mat/020814. 

Appendix: Appendix 

This appendix summarizes, without derivations, useful 
technical results that arc alluded to in the main text. 



A.l Kubo formula in terms of path integrals 

The Kubo formula for the DC conductivity of a d- 
dimensional conductor can be expressed as 



ctdc = - lim - V / 
w ^o d ' J 



dx 2 in' 



22 



(<"o) 



G 



11', 22 



'(W ) 



J — ( 



dt 12 e l ^9{t 12 ) c {1 



(AT) 



C[ll',22'] = ^{[$\h, X V )$(tl, Xl),ft (t 2 , Xytyfa, X 2 )]) H , 

where j\y = (^i ~ ^l') an d a uniform applied 
electric field E{lo q ) was represented using a scalar po- 
tential, V cxt {uJo,r 2 ) = —r 2 ■ E(ll> )- A path integral 
representation for Cjn' 22'] can be derived using the fol- 
lowing strategy, adapted 18 ^ from GZ's Rcf.: 2 ' (1) intro- 
duce a source term, in which an artificial source field 
«2'2 couples to , 4>^(t 2 ,x 2 i)'ip(t 2 ,x 2 ), and write C , [n' ! 22 / ] as 
the linear response of the single-particle density matrix 
piv — (V^(*i) xv)^){t\, xx))h to the source field V22' ■ (2) 
Decouple the interaction using a Hubbard-Stratonovitch 
transformation, thereby introducing a functional integral 
over real scalar fields V F / B , the so-called "interaction 
fields" ; these then constitute a dynamic, dissipative en- 
vironment with which the electrons interact. (3) Derive 
an equation of motion for p\ v , the single-particle den- 
sity matrix for a given, fixed configuration of the fields 
Vp/R, and linearize it in -52'2, to obtain an equation of 
motion for the linear response 5p\ v {f) to the source field. 
(4) Formally integrate this equation of motion by intro- 
ducing a path integral J V{RP) over the coordinates 
and momenta of the single degree of freedom associated 
with the single-particle density matrix SpYy ■ (5) Use the 
RPA-approximation to bring the effective action Sy that 
governs the dynamics of the fields V F / B into a quadratic 
form. (6) Neglect the effect of the interaction on the 
single-particle density matrix wherever it still occurs, i.e. 
replace p£ by the free single-particle density matrix 

p% = (^(xS( Xi )) =J2rx(xj)M^)no(tx) , (A-2) 



where (6)0 = Tre[-^ O]/Tr[e^ ]. (7) Perform the 
functional integral (which steps (5) and (6) have ren- 
dered Gaussian) over the fields V F / B ; the environment is 
thereby integrated out, and its effects on the dynamics of 
the single particle are encoded in an influence functional 
of the form e~( lSR+Sl \ The final result of this strategy 
is that C[n/.22'] can be written as 17 ) [cf. (11.49)] 



lim 

to^ — 00 



/ 



dx 2 {^l2',22,21'(^l)^2;io) — ^12,22',21'(il)*2;io)}) 



'12', 22', 21' 



{tl,t2;to))y 



V(RP) 



(A-3) 



•2f,2' h 



Of, b 



V{RP)- / dx dxo fa e -liSR+Si](ti,t )/h 
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where = i.e. the indices 2 and 2 refer to the 
same time £2, but the integration variable £2 is indepen- 
dent of X2- The symbol J V(RP) stands for the double 
path integral defined in Eq. (4). The complex weight- 
ing functional e l ^° •* occuring therein involves the 
action for a single, free electron, in the mixed coordinate- 
momentum representation, 

S^t^R^t^P^t,)} (A-4) 

= t dt 3 [P a (t 3 ) ■ d t3 R a (t 3 ) - h a (R a (t 3 ),P a (t 3 ))] , 
Jti 



ho(R a ,P a ) 



2m 



+ V imp (R a ) - p • 



(A-5) 



The first line of Eq. (A-3) corresponds to Eq. (3), and 
its significance as describing Cooperon propagation along 
two paths from ry at time ti to T\ at t\, and n> at time 
ti to Ti at t2 1 is explained after Eq. (4). The integrals in 
the second line of Eq. (A-3) "initialize" this Cooperon 19 -* 
at an initial time to ~~ > —00, by including forward and 
backward propagation [see the left (non-loop) part of 
Fig. 1] between the times to and ti, from and to some 
initial positions r and rg, which are weighted by the ini- 
tial, free single-particle density matrix p ^ of Eq. (A-2). 

The effective action iSjj+Sj in the influence functional 



-(zSr+Si) 



depends via Eq. (6a) on the single-particle density matrix 
in a mixed coordinate-momentum representation, 20 -* 

p F {R F , P F ) = J dr F e -^ F - 5rF p°(R F + 5r F , R F ) , 

(R B , P B ) = J dr B e -^ B - SrB p°(R B ,R B - 5r B ) 

(A-6) 



Po 



where p {ri,rj) — p^ is defined in Eq. (A-2). 

Furthermore, iSr + Si depends on two purely real 
functions, R(t,r) and I(t,r), that are defined as follows 
via their Fourier transforms [cf. (11.56,57)]: 



(R/I)(t,R) = 
R(u,k) = 



dkduo 
V int (k) 



-iojt-\-ik-R 



(R/I)(co,k) , 



(A-7a) 
(A-7b) 



1- V int (fc)x(w,fc) 
I(uj, k) = - coth(fiw/2T) lmR(oj, k) . 

Here V lnt (q) and x(o>,fc) are the Fourier transforms of 
the interaction potential V" mt (|ri — r 2 |) and the charge 
susceptibility X12, which can be written as 



X12 



-ie 2 ft(Gf 2 G£+G&G&) . (A- 



A. 2 Derivation of Cooperon self energy from influence 
functional 

The path integrals / V(RP) of Eq. (3) can be given 
a precise definition in terms of the standard time-slicing 
procedure for path integrals, with one coordinate and 
one momentum integral for every time slice. 20 ) For each 
time slice (labeled by n, say), the momentum integral 
/ dP n can then easily be performed, since it simply has 
the effect of converting the expressions occuring in the 
action for that time slice from the mixed coordinate- 
momentum representation to the coordinate-only repre- 
sentation. Thus, the free action Sft[R,P] is mapped to 

§S(t,t')[R a ,R a ] = £dt 3 [\mR a2 {t 3 ) - V imp (R a (t 3 )) 

the density matrix at time slice n is converted 20 - 1 from 
pg(K,P„ a ) to /5°(r», »•»_!), and C 12 ,, 21 , of Eq. (3) can 
be rewritten as 



12', 21' 



f V'(R) e-liSn+SiKtutoyn ; ( A . 1Q ) 

J 2 f ,2b 



where the integral / V\R) is used as a shorthand for 



V\R) 



in Eq. (A-3) is given by Eqs. (5) and (6). S R J2 'f^b 



R F (t!)=n 



R"(t 2 )=r 2 , 

H B (tl)=TV 



V'R* (t 3 ) 



(ATI) 



x / V'R B (t 3 )e^~ s "^ h 

lR B (t 2 )=r 2 



The effective action iSr + Si in Eq. (A- 10) is found to 
have the following form 21 ' [with Su = 8 a - ai 5{ri — 7^) and 
{R/I)- la j a , = (R/I)(U - tj,r?(u) - rf(tj))}: 



SR/i(ti,t 2 ) = J2 [ dt 3a f 
Jt 2 Jt- 



1 dt 3a [ dt , .!.'< I , , (A-12) 
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G^ A ^ K are the retarded, advanced and Keldysh Green's 
functions of the nonintcracting, disordered system: 

Gf/ A = T \0(±U 3 )Y j r x {x 3 )M^y- lixUl,h . (A-9a) 

A 

= -^^M^e-^^l 2n (6)] - 



Now expand the effective action in powers of S R /i, 

If.I'b 



Cl2',2l' = E W /, - V 'W 
N=0 ' ^2 F ,2 B 



(A-14) 
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and analyze the structure of this expansion, which can 
be viewed as a type of Dyson equation for the Cooperon. 
The N = term simply yields free propagation between 
the specified end points, C° 2 , 2l , = K 2 G r > 1f2 'fG£ bV . 
The N = 1 term of Eq. (A- 14) turns out to have the 
form of a forward-backward propagator sandwidching 
self-energy insertions, 21 ^ say that are depicted di- 

agrammatically in Fig. 2 and are given by: 

f lF,1 ' B ~RII f 

/ V\R)Lj a , (tzM)= I dx3 F dx3 B dx4 F dx iB 
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as p" 
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Eoo 
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[-av ri )'/2m+y imp (r,)- 



2T 



Now, GZ neglect the fact that SV ri does not commute with 
Vi m p(ri), arguing that this is justified here in the semiclassi- 
cal limit h — > 0. We disagree, since in the above expression for 
p% every neglected power of h is accompanied by at least one 
power of 1/T, i.e. the neglected terms are of order hj (Tr^^), 
where T c h ar is a time scale characteristic of the disorder po- 
tential (rough estimates suggest r c h a r — r e i). Thus, they can 



be neglected only for T > h/r e i . 
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n^nf^ 1 (/<*<) llf =1 (/ J$*y<~ ) , «*c. In doing 
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If T > 5J_ 1 and r™ are the end-of-interval coordinates for time 
slice n, and <5r™ = — is the relative coordinate [which 

is integrated over in the Fourier transforms of Eq. (A-6)], then 
we choose the "center-of-mass" coordinate as R£ = r^_j and 
= , respectively. 

21) For each occurrence in Eq. (AT3) for L R//I of a pair of indices, 
one without bar, one with, e.g. 4 a and 4 a , the corresponding 
coordinates rj and r| are both associated with the same 
time <4, and integrated over, f dr^dr^, in the path integral 

/ V\R). When evaluating Eq. (A-15) for L 1 or L R , one of 
these integrals can be used, respectively, to collapse the 5(rJ — 
function in L 1 , or to construct a Kcldysh function from 
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tional if their arguments are real, u i§u can never cancel any 
contribution from Sj" [discussion before (III. 22)]. However, 
this argument overlooks the fact that Sr and Si are both 
functionals of the path integral variables R(t3) and P(t:j), 
and that the integration measure e l (^o ) in J T>(RP) is 
complex (GZ's argument would be correct if the integration 
measure were real). 

23) However, the analysis there (critized in cond- mat/0208264) is 
yet again semiclassical, and hence cannot properly deal with 
the intrinsic problem (IP) pointed out in section 5, namely 
that the density matrix becomes a sharp function at low tem- 
peratures. Thus, we believe that the expansions in Eqs. (40) 
to (42) of cond-mat/020814 will become singular at low T. 



